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ABSTRACT 

Least  squares  methods  are  used  with  NIMBUS  IV  SIRS-B 
clear -column  radiances  as  independent  variables  and  two 
identical-period  sources  of  geopotential  heights  (at  the 
500  mb  and  300  mb  pressure  levels)  as  dependent  variables. 
Regression  equations  are  developed  for  each  of  the  sets  of 
geopotential  heights  for  three  latitude  bands  in  the 
Northern  Hemisphere.   Three-day  and  four -day  data  bases  were 
used  on  each  regression  specification.   Each  set  of  equa- 
tions was  tested  on  independent  data  at  a  subsequent 
composite  period  of  12-hours  and  24-hours  following  the 
dependent  pooled  samples.   The  regression  results  of  both 
the  dependent  and  independent  tests  were  examined  for 
possible  operational  usefulness. 
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I.   INTRODUCTION 

The  Satellite  Infrared  Spectrometer  (SIRS-B)  carried  on 
the  NIMBUS  IV  spacecraft  system  relays  to  National 
Aeronautics  and  Space  Administration  (NASA) ,  at  Goddard 
Space  Flight  Center,  radiance  observations  taken  in  seven 
distinct  channels  of  the  15  |jm  band  of  carbon  dioxide,  and 
one  at  899.0  cm   in  the  atmospheric  window.   In  addition, 
there  are  six  channels  in  the  water-vapor  rotational  band 
also  being  relayed  by  NIMBUS  IV,  but  these  latter  six 
channels  will  not  bear  on  the  subject  of  this  dissertation. 

Each  of  the  eight  primary  channels  (699.3  cm-1, 
677.8  cm   ,  692.3  cm   ,  699.3  cm"1,  706.3  cm"1,  714.3  cm"1, 
750.0  cm   ,  and  899.3  cm   )  is  sensitive  to  temperature 
variations  in  the  tropospheric-stratospheric  range  of  heights 
In  addition,  each  channel  has  a  maximum  transmissivity  in  a 
discrete  height  range  above  the  surface.   Smith,  Woolf,  and 
Jacob  (1970)  presented  a  graph  depicting   d T  /dlnp   versus 
pressure  where   t   is  the  transmissivity  and  the  derivative, 
dT/dlnp,  becomes  a  weight  factor  of  the  relative  contribu- 
tions of  the  various  elevations  to  the  total  emitted 
radiance  in  the  channel.   It  can  be  seen  in  Figure  1  and 
Equation  (1) ,  after  Smith  et  al  (1970) ,  that  the  different 
channel  radiances  are  most  responsive  to  changes  in  tempera- 
ture which  occur  at  the  level  of  their  peaks  in  dT  /din  p. 
Each  channel  tends  to  signal  the  occurrence  of  temperature 


12 


changes  at  or  near  its  own  pressure-level  peak,  as 
discernible  from  Figure  1. 


0.6     0.8 
dT/dlnp 


Fig.  1  Derivative  of  transmittance  with  respect 
to  the  logarithm  of  pressure  for  eight 
SIRS-A  channels  [after  Smith  et  al 
(1970)]  . 

The  radiative  transfer  integral  applied  to  a  clear 
column  of  air  for  which  the  sounding   T(p)   is  known  gives 
the  radiance  as 


N 


/(PS) 
B  [vi.T(p)]  dTi 


(1) 


where   B  \/i/T(p)    is  the  Planck  radiance  at  wave  number 
VV*  and   T(p),  the  temperature  at  pressure  level   p.   In 
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Equation  (1)  the  subscript  "s"  indicates  a  surface  value  and 
T(Vi'  P)  is  tne  fractional  transmittance  of  the  atmosphere 
in  the  spectral  interval  centered  at   Vj^   from  pressure 
level  p   up  to  the  satellite.   If  the  sounding  is  known, 
the  radiance  in  the  channel  centered  at  wave  number  y/. 
can  be  computed  since   t(  y/\,    p)  is  known  as  a  function  of 
pressure  p    Smith  (1970)  .   The  inverse  problem  has  become 
more  relevant  since  the  advent  of  the  SIRS  instrument, 
namely  in  the  use  of  observations  in  the  eight  channels  to 
deduce  the  clear  atmosphere  structure   T  =  T(p) .   The  usual 
method  of  solving  an  integral  equation  for  the  function 
B(  Vi'  T)  was  found  to  be  unstable,  but  it  was  found    Smith 
et  al  (1970)   that  regression  methods  relating  temperatures 
and  geopotential  heights  at  standard  pressure  levels  were 
quite  well  described  by  the  results  of  SIRS  scan-soundings 

Smith,  Woolf,  and  Jacob  (1970)  found  it  convenient  to 
convert  the  clear  column  radiances  into  equivalent 
"brightness  temperatures,"   Tg (  y/j) .   The  definition  of 
TB  (  y/±)       i-s  based  upon  the  particular  Planck  function  for 
temperature   TB (  y/^)      which  equals  the  clear-column  channel 
radiance  N(  V^) ,  that  is 


Wv  tb  <Vi>]  =  »cv±> 


Byr  tb  (Vi)      = 

so    that  by   solving    for      Tg     we  have 


Tb(V±)    =   C2Vi     |  In    J    C1Vi3/N(>/i)l     +    1      > 


-1 

(3) 
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Here  the  constants   C^   and   C2  have  the  values 

C±  =   1.19061  x  10   erg-cm  -sec  x  -steradian 

C2  =  1.43868   cin-0K 
Smith  et  al   (1970)  derived  regression  equations  for  contour 
heights  at  standard  pressure  levels  (see  Table  1)  utilizing 
the  eight  simultaneous  brightness  temperatures   TB(\Z^)   as 
independent  variables.   As  dependent  variables  they  employed 
a  large  sample  of  radiosonde-reported  geopotential-heights 
at  the  mandatory  levels.  "Their  results  for  the  month  of 
September,  1969,  using  a  two-week  sample,  was  stated  in  the 

nonlinear  regression  form 

8 
Z(Pj)  =  Z(Pj)  + '  ]^b  (vVPj).[TB(v'i)  "  VVj.)]  + 


i  =1 
8 


(4) 


£b,(Vi.Pj>  [Vi)  "  *B<V±>] 


i  =  1 

with  an  analogous  regression  equation  for   T(p^) .   Their 
regression  coefficients  for  the  September,  1969  SIRS-A 
example  are  shown  in  Table  1,  where  the  coefficients  for 
T(pj)  are  also  listed.   For  the  data  period,  24-29  December 
1970,  employed  in  this  study,  an  analogous  set  of  regression 
equations  based  upon  updated  SIRS-climatology  would  be  used. 

It  is  clear  that  the  form  in  Equation  (4) ,  as  employed 
by  Smith  et  al  (1970) ,  essentially  builds  from  the  tempera- 
ture structure  of  the  remotely-sensed  atmospheric  columns 
viith  ps  =  1000  mb  and  reconstructs  a  height  profile  by 
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statistically  relating   Z(pj)   through  the  hydrostatic 
equation  to   Z  (p-^)   at  any  other  level  desired. 

The  objective  of  this  study  was  to  obtain  statistical 
relationships  for  geopotential  heights  of  the  500  mb  and 
300  mb  pressure  levels  by  stepwise  regression  analysis 
between  these  variables  and  the  eight  clear-column  satellite 
radiances.   Further,  it  was  desired  to  do  this  without  prior 
computation  of  the  brightness  temperatures.   The  latter 
procedure  is  an  unnecessary  step  if  corrected  radiances  are 
known,  as  was  the  case  in  this  study.   If  a  sixteen- 
predictor  equation  were  developed,  Equation  (4)  would  then 
be  modified  to  include  the   Nj_   and  Nj_    terms  (i  =  1,...8) 
with  appropriate  coefficients. 

With  the  use  of  NIMBUS  IV  SIRS-B  corrected  radiances, 
Ni,  and  National  Environmental  Satellite  Center  (NESC) 
regression-computed  contour  heights  at  the  same  geographic 
locations,  diagnostic  regression  equations  for  the  500  mb 
and  300  mb  geopotential  heights  were  developed.   Separate 
regression  equations  were  developed  using  both  three-day 
and  four-day  pooled  data  bases.   The  regression  equations 
developed  were  tested  for  validity  on  an  independent  sample 
in  each  case. 

The  same  procedure  was  applied  to  Fleet  Numerical 
Weather  Central  (FNWC)  analyzed  geopotential  heights  at  the 
500  mb  and  300  mb  levels  for  the  same  positions  and  times. 
In  the  latter  case,  the  contour  heights  had  in  no  way  been 
related  to  SIRS-B  scan-soundings.   Thus,  any  diagnostic  and 
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prognostic  relationships  discovered  here  would  have  to  be 
based  upon  the  statistical  correlation  coefficients 
R  [z(5/3),Ni]  '  i  =  1#...#8/  as  previously  suggested  by 
Smith  et  al. 
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II.   DATA  PROCESSING 

A.   THE  ORIGINAL  NESC  DATA 

The  original  data  obtained  from  the  National  Environmental 
Satellite  Center  (NESC)  consisted  of  a  computer  listing,  by 
days,  of  selected  SIRS-B  scan  spots  during  the  period 
Dec.  1970  -  Jan.  1971.   Included  in  this  listing  were  the 
geographical  sites  and  times  of  each  observation,  amounts 
and  heights  of  effective  cloud  cover,  temperatures  and 
heights  of  the  standard  levels,  and  a  double  set  of 
radiances.   The  latter  were  the  original  set  of  uncorrected 
radiances  and  the  cloud-corrected  values  for  each  channel. 
The  uncorrected  radiances  in  the  window  channel  were  missing 
since  that  channel  had  become  inoperative  by  December,  1970. 
The  sample  of  scan-spots  presented  by  NESC  was  therefore 
restricted  to  sites  having  known  surface  temperatures  so 
that  a  realistic  effective  window  channel  radiance,  B(v"/Ts) 
for  a  surface  pressure  level  ps  =  1000  mb  could  be  assigned 
for  each  location. 

It  should  be  pointed  out  at  this  stage  that  the  SIRS-B 
radiometer  carried  aboard  the  NIMBUS  IV  differed  from  the 
SIRS-A  radiometer  aboard  the  NIMBUS  III  in  that  the  SIRS-B 
radiometer  had  a  "side-looking"  capability  as  well  as  the 
vertical  mode.  This  permitted  a  wider  range  of  soundings 
in  SIRS-B  than  was  possible  with  SIRS-A.  These  side-view 
radiances  were  normalized  to  a  single  atmospheric 
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emitter-depth  located  over  the  scan  spot  so  that  all  the 
data  was  compatible,  apart  from  variations  in  cloud  cover 
and  surface-layer  lapse  variations. 

Because  of  the  scan-spot  area  of  (225  km)  ,  most  spots 
viewed  were  partially  cloud-filled  with  some  effective 
fraction,  A.   This  fact  complicated  the  determination  of  the 
cloud-corrected   T(p)   profiles  since  the  radiances  were 
reduced  to 


N(V±)    -  A   J    B       jVi<VPc>]      t(VVPc) 
"    J        B     [  Vi/T(p)  ]     d  Ti| 


1 
+ 


(5) 


(1    -A)    J     BJV^Tte^J     T  (VVPS) 

J  B         [VrT(p)]    d   tJ 

assuming  a  one-level  cloud  model  with  the  top  at  pc.   Here 
the  part  of  Equation  (5)  within  the  second  brace  is  that 
part  of  NCv'-j)   which  represents  the  desired  uncontaminated 
radiance  in  channel   i  which  escapes  through  the  cloud-free 
fraction,  1  -  A. 

Smith  et  al  (1970)  have  outlined  a  method  for  determining 
a  cloud-contamination  correction  which  yields  the  estimate 
for  corrected  radiance   Nc(\/^)   in  each  channel.   Since  the 
work  in  this  study  was  dependent  upon  the  accuracy  of  the 
corrected  radiances,  a  brief  explanation  of  this  correction 
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procedure    is   given,    based   on   the   one-level   cloud   model  with 
top   at   level      pc.       Equation    (5)    can  be  broken    into   two  parts, 

Nc(Vi)     =   N(  V±)    +    CfvA)  (6) 

where  the  correction   C(V^)  is  designed  to  eliminate  the 
effects  of  cloud-reduction  of  the  radiance,  and  may  be 
obtained  by  dividing  Equation  (5)  by  the  factor  1-A.   How- 
ever, this  will  not  be  a  permissible  procedure  since  the 
cloud  cover,  A,  may  be  very  close  to  unity.   Instead,  Smith 
et  al  (1970)  suggest  an  iterative  procedure  whereby   C (  \/ . ) 
of  Equation  (6)  may  be  obtained  in  the  form 


cCv^)  =  a  |  Nc(Vi)  -  x|yrPc,T(P)]  J 


(7) 


The  function  X  V-»P  'T(p)    represents  the  contents  of  the 
second  brace  of  Equation  (5) .   Successive  corrections 
C(V.)   are  then  determined  for  each  channel  4,...., 8,  by 
making  the  assumption  that  the  radiances  observed  in 
channels  2  and  3  (7  50.0  cm   and  714.0  cm"  )  are  correct 
during  iterations  while  correcting  the  remaining  set  of 
radiances  N^,...Ng. 

In  determining  a  first  guess  to  a  clear  column  atmosphere 
T(p)  the  five  most  opaque  channels  in  Figure  1  are  assumed 
correct  at  the  first  iteration,  as  well  as  the  pre- 
determined window  channel.   The  NESC  correction  scheme  then 
makes  use  of  a  six-predictor  regression  equation  for   T(p) 
(analogous  to  Equation  (4)  but  lacking  channel  2  and  3  terms) 
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which  then  specifies  the  required  first  guess   T(p-j)   at  all 
standard  levels.   With  this  profile  now  specified  the 
best-fit  consistent  choices  of  A,  1-A,  and  p   which  comes 
closest  to  providing  the  best-fit  solutions   N(  y/^)       to 
Equation  (5)  for  channels  2  and  3  are  obtained.   In  this 
test  pc   is  subject  to  variation  over  the  standard  levels 

p  =  200,  250, 850  mb   and  the  best-fit  A,  1-A,  are 

thereby  chosen.   At  each  iteration,  there  is  an  improved  set 
of  six  corrected  radiances,  which  leads  to  an  improved 
specification  T(p^).   This  then  leads  to  improved  newly 
corrected  radiances   Nc (  s/-)     (i  =  4,....,8)  in  accord  with 
Equations  (6)  and  (7) .   Finally,  in  a  limited  number  of 
iterations   C(v'.j)   approaches  a  convergent  value.   At  this 
final  stage,  the  six-predictor  regression  equation  for 
T(Pj)  gives  estimates  of  corrected  radiances  in  channels  2 
and  3.   The  listed  statistics  relating  pc   and  A  are  of 
no  particular  importance  since  these  quantities  relate  only 
to  an  effective  cloud  cover  over  an  area  (225  km) 2,  and  the 
cloud-depth  is  considered  irrelevant  in  obtaining   N^y'^). 
The  important  factor  to  be  tested  in  this  paper  is  the 
concept  that  geopotential  heights  as  taken  from  a  large 
radiosonde  sample  are  highly  correlated  with  the  set  of 
clear-column  radiances,  i.e.,  to  the  estimated   T(p) 
structure  of  the  atmosphere.   Thus  the  NESC  SIRS-B  data 
obtained  by  the  December  1970  version  of  the  NESC  regression 
equation  system  provided  also  the  geopotential  heights  and 


22 


temperatures  at  all  standard  levels  in  the  atmosphere  from 
p  =  1000  mb  to  10  mb. 

In  developing  the  NESC  regression  equations,  Smith  et  al 
(1970)  utilized  a  data  sample  of  approximately  two  weeks, 
continually  updating  the  sample  period  closer  to  the  time 
when  the  independent  calculations  were  made.   The  heights 
used  in  their  regression  equations  were  interpolated  in 
space  and  time  from  National  Meteorological  Center  analyzed 
fields  at  standard  pressure  levels  to  the  geographical  loca- 
tions and  observed  times  of  the  SIRS  measurements.   Only 
portions  of  analyses  and  NIMBUS  III  data  close  to  radiosonde 
observations  were  used  for  developing  their  dependent  data 
regression  coefficients  giving  heights  from  corrected 
radiances.   These  NESC  dependent  samples  were  divided  into 
latitude  bands  (18-40°N,  35-55°N,  and  50-80°N) .   A  similar 
technique  was  used  by  NESC  in  making  available  the  data  from 
December  1970,  a  portion  of  which  was  used  in  this  thesis. 
This  data  consisted  of  corrected  radiances  at  scan-spots 
which  were  used  in  conjunction  with  an  eight-predictor 
regression  equation  for  generating  temperature  and  height 
values  for  the  standard  levels  from  the  SIRS  observations. 
These  computed  heights  and  temperatures  appeared  on  the  NESC 
data  listings  for  24-28  December  1970,  which  were  used  in 
this  study. 

The  NESC  regression-derived  heights  of  the  500  mb  and 
300  mb  levels,  position,  time,  and  clear-column  radiances 
for  each  sounding  from  24-27  December,  1970,  were  used  for 
-the  regression  analysis  described  in  Chapter  III. 
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B.   FNWC  DATA 

Height  fields  for  0000  GMT  and  1200  GMT  covering  the 
period  24-29  December  1970,  were  obtained  from  the  data- 
tapes  of  Fleet  Numerical  Weather  Central  (FNWC)  at  Monterey 
for  both  the  500  mb  and  300  mb  levels.   The  contour  heights 
were  interpolated  in  space  and  time  to  the  same  positions 
and  times  of  the  satellite  data  obtained  from  NESC.   The 
spatial  values  of  geopotential  height  were  computed  by  use 
of  a  double-quadratic  Bessel  interpolation  scheme  applied  to 
both  the  synoptic  times  preceding  and  following  each  satel- 
lite observation.   Then  the  heights  were  interpolated 
linearly  to  the  time  of  the  satellite  observation.   The 
computer  program  used  to  accomplish  the  interpolation  scheme 
is  listed  in  the  Appendix. 

The  FNWC  data  provided  an  independent  set  of  geopotential 
heights  for  the  500  mb  and  300  mb  levels  from  which  a  second 
set  of  regression  relations  could  be  derived  for  comparison 
with  the  regression  analyses  performed  on  the  NESC  data. 
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III.   DEPENDENT  DATA  REGRESSIONS 


A.   DEVELOPMENT  OF  REGRESSION  EQUATIONS  USING  NESC  GENERATED 
VALUES  OF  500  MB  AND  300  MB  GEOPOTENTIAL  HEIGHTS. 
The  NESC-generated  300  mb  and  500  nib  pressure  levels 
(ZN(5)   and   znO)   respectively)   were  selected  for  investi- 
gation with  the  corresponding  500  mb  and  300  mb  FNWC 
analyses,  which  presently  have  no  SIRS  input.   A  multiple 
linear  regression  analysis  technique  was  used  in  this  study 
corresponding  to  the  8-predictor  scheme  of  Smith,  Woolf,  and 
Jacob  (1970),  but  containing  no  nonlinear  terms.   The  speci- 
fic program  used  was  BIMED  02R,  a  version  of  the  stepwise 
linear  regression  technique  currently  on  file  at  the 
Statistical  Library  of  the  W.  R.  Church  Computer  Facility 
at  the  Naval  Postgraduate  School,  after  Dixon  (1966) .   The 
result  of  this  technique  when  applied  to  simultaneous  data 
samples  of  the  form  [Z^N-j^N^,  .  .  ,Ng]   is  to  produce 
regression  equations  of  the  form 


8 


=  Co  +  X)  Ci 


N-;  -   (C0,  C-.  ,....,  Co  ) 


i=l 


(8) 


Here  the  subscripts  1,..,8  indicate  the  channel  numbering 
system  and  are  also  used  for  variable  numbering.   A  series 
of  dependent  test  regressions  were  performed  on  both   Zjj(5) 
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and   ZN(3)  using  a  pooled  three-day  data  base  and  a  pooled 
four -day  data  base  of  simultaneous  values  of  SIRS-B  correc- 
ted radiances.   Due  to  the  nature  of  the  data  available, 
each  data-day  presented  a  composite  of  scan-soundings 
centered  about  1200  GMT,  varying  up  to  six  hours  on  either 
side  of  this  central  time.   The  data  was  further  divided 
into  three  overlapping  latitude  bands.   The  limits  of  each 
band  are  listed  below: 

Band  1.    20°N  -  40°N 

Band  2     35°N  -  55°N 


Band  3     50°N  -  70°N 


The  resulting  coefficients  (CQ,  C1,....,Cg)  for  the  three- 
day  and  four -day  regression  equations  are  listed  in  Tables  2 
and  3  for  the  500  mb  and  300  mb  dependent  tests, 
respectively. 
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TABLE  2.   500  mb  NESC  Dependent  Regression 

Coefficients  for  Three-  and  Four-Day 
Data  Base,  Giving  Best-fit  for  Eq.  (8) 


3-Day  Sample  Base 


Band  1 

Band  2 

Band  3 

Sample  Size 

171 

190 

177 

Co 

5473.77 

5689.80 

4620.48 

cl 

6.66 

-2.80 

-2.90 

c2 

-17.56 

-3.32 

-11.68 

c3 

23.85 

41.83 

62.59 

c4 

28.  Q6 

-10.83 

-17 .  87 

c5 

-62.22 

-66.78 

-36.90 

C6 

-6.80 

4.02 

-4.14 

C7 

21.20 

19.62 

11.74 

C0 

-0.33 

-1.40 

* 

4-Day  Sample  Base 


Sample  Size 


C3 


-8 


Band  1 

Band  2 

Band  3 

171 

227 

216 

5497.55 

5756.70 

4555.32 

6.53 

-1.18 

-4.45 

-15.98 

-7.7  5 

-6.38 

20,13 

45.24 

54.55 

34.33 

-12.75 

-10.93 

-67.57 

-64.77 

-35.75 

-5.88 

2.64 

-7.94 

24.18 

19.54 

1.00 

-3.26 

-1.13 

1.45 

Indicates  an  insignificant  contribution  to  the  analysis 
of  variance  by  this  predictor. 
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TABLE  3.   300  mb  NESC  Dependent  Regression  Coefficients 
for  Three- and  Four -Day  Data  Base,  Giving 
Best  Fit  for  Eq.  (8) 


3-Day  Sample  Base 


Band  1 

Band  2 

Band  3 

Sample  Size 

171 

190 

177 

co 

8976.59 

8490.63 

7166.52 

cl 

12.79 

4.00 

1.57 

c2 

-22.36 

-6.08 

-15.54 

C3 

4.04 

21.62 

54.11 

c4 

87.54 

47.07 

23.71 

C5 

-105.20 

-96.10 

-48.47 

C6 

-15.77 

4.64 

-4.34 

C7 

23.90 

15.88 

6.14 

c8 

5.67 

1.62 

1.29 

4-Day  Sample  Base 


Band  1 

Band  2 

Band  3 

Sample  Size 

206 

227 

216 

Co 

8986.32 

8613.84 

7102.40 

Cl 

13.26 

6.65 

* 

C2 

-21.76 

-13.75 

-10.18 

C3 

* 

29.09 

45.98 

c4 

95.49 

41.53 

30.92 

C5 

-110.81 

-93.15 

-47.49 

c6 

-14.34 

3.13 

-8.48 

C7 

28.03 

16.26 

5.37 

C8 

1.25 

1.65 

2.94 

* 

Indicates  an  insignificant  contribution  to  the  analysis 

of  variance  by  this  predictor. 


28 


The  coefficients  in  Tables  2  and  3  are  considerably 
different  in  each  latitude  band  which  tends  to  justify  the 
use  of  separate  latitude  bands  in  the  analyses.   Whether  the 
limits  of  each  latitude  band  are  optimum  or  not  is  undeter- 
mined, but  it  was  felt  that  any  further  reduction  in  the 
width  of  each  band  would  reduce  the  number  of  soundings  used 
in  each  regression  equation  and  thereby  have  a  detrimental 
effect  on  a  statistical  analysis.   The  differences  between 
the  corresponding  coefficients  from  the  three-  and  four -day 
sample  regressions  were  generally  insignificant;  however, 
this  conclusion  will  be  further  tested  by  examination  of  the 
fractional  explainted  variance  accounted  for  by  the  appro- 
priate regression  equations  when  applied  to  the  data- 
samples. 

In  both  Tables  2  and  3,  the  coefficients  corresponding 
to  channels  3,  4,  and  5  appear  to  be  weighted  more  heavily 
than  any  other  triad  of  coefficients,  which  appears  to  be 
consistent  with  the  response  levels  of  Figure  1. 

A  useful  statistic  for  the  fractional  unexplained 
variance  is 

2 

o        n  -  k  -  1     /std.  error  \ 
1  "  r2  -  -T-T-i—       ^5   ho„  (9) 


n  -  k  -  1     / std.  error  \ 
n  -  1       I  std.  dev.  J 


where   k   is  the  number  of  predictors  and   n   the  sample 
size. 
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Thus  the  fractional  explained  variance  is  easily 
obtained  from  Equation  (9)  as  v 

r2  m        _   n  -  k  -  1    /std   error  \2 

n  -  1      I   std.  dev.  J 

Values  for   R^   and  of  standard  errors  for  each  of  the 
dependent  tests  corresponding  to  Tables  2  and  3  are  listed 
in  Table  4. 

TABLE  4.   R2  and  Standard  Errors  (gpm)  for  NESC 
Dependent  Regression  Equations  for 
500  mb  and- 300  mb  Three-  and  Four -day 
samples. 


3 -Day  Sample      Band  1  Band  2         Band  3 


500 

mb 

Band  1 

Band  2 

.8050 

.8222 

61.72 

77.13 

R2  .8050  .8222  •         .9050 

Std.  error (m)      61.72  77.13  69.68 


4 -Day  Sample 

R2  .8039  .8154  .9051 

Std.  error  (m)      60.55  76.78   %.        69.79 

300  mb 

3 -Day  Sample      Band  1  Band  2         Band  3 

R2  .8790  .8693  .9351 

Std.  error(m)      77.51  91.80.  75.36 

•\  . 
4 -Day  Sample 

R2  .8764  .8622  .9349 

Std.  error(m)      76.22  91.42  75.50 
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As  indicated  in  Table  4,  Band  3  appears  to  produce  higher 
values  for   R  ,  and  therefore  better  results,  than  the  other 
two  bands  in  the  NESC  dependent  regression  analysis.   The 
standard  errors  for  the  300  mb  tests  were  generally  higher 
than  the  500  mb  tests,  but  one  might  expect  this  to  happen 
since  the  standard  deviations  of  zi$(3)    are  greater  than 
those  of   ZN(5).   In  general,  the  300  mb  tests  gave  higher 
values  of  R2   upon  specification  than  the  corresponding 
500  mb  tests. 

B.   DEVELOPMENT  OF  REGRESSION  EQUATIONS  USING 

FNWC- INTERPOLATED  VALUES  OF  500  MB  AND  300  MB 
GEOPOTENTIAL  HEIGHTS. 

Procedures  similar  to  those  used  in  Part  A  of  this 
chapter  were  used  on  the  FNWC  500  mb  and  300  mb  geopotential 
heights  *ZF(5)   and   Zp(3)   respectively]  ,  corresponding 
to  the  identical  SIRS  scan-sounding  positions  and  times. 
The  same  latitude  band  limits  as  specified  in  III. A.  were 
used  for  both  the  three-  and  four -day  regression  analyses. 
The  resulting  coefficients  for  the  FNWC  regression  equations 
are  listed  in  Tables  5  and  6  for  the  500  mb  and  300  mb 
dependent  tests,  respectively. 
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TABLE  5.   500  mb  FNWC  Dependent  Regression  Coefficients 
for  Three-  and  Four -day  Data  Base 

3 -DAY  SAMPLE  BASE 


Band  1 

Band  2 

Band  3 

Sample  Size 

171 

190 

177 

co 

5617.09 

5583.51 

4615.09 

cl 

7.58 

-2.25 

-2.50 

c2 

-19.82 

-2.58 

-9.43 

c3 

23.50 

38.35 

57.19 

C4 

25.55 

-6.14 

-16.06 

C5 

-53.38 

-67.11 

-32.17 

c6 

-13.71 

5.85 

-9.10 

c7 

23.14 

19.19 

12.22 

C8 

-1.35 

-2.37 

0.80 

4 -DAY  SAMPLE  BASE 


Band  1 

Band  2 

Band  3 

Sample  Size 

206 

227 

216 

Co 

5653.75 

5762.27 

4595.81 

Cl 

7.70 

0.36 

-2.24 

C2 

-20.43 

-10.14 

-9.16 

C3 

24.21 

46.38 

55.72 

c4 

28.86 

-13.41 

-15.03 

C5 

-59.44 

-64.52 

-30.50 

c6 

-12.29 

4.74 

-9.58 

c7 

27.24 

19.12 

8.71 

C8 

-5.01 

-2.15 

3.26 
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TABLE  6.   300  nib  FNWC  Dependent  Regression  Coefficients 
for  Three-  and  Four-day  Data  Base 

3 -DAY  SAMPLE  BASE 


Band  1 

Band  2 

Band  3 

Sample  Size 

171 

190 

177 

Co 

9156.82 

8451.37 

7337.29 

Cl 

14.03 

2.99 

* 

C2 

-33.45 

-8.12 

-18.36 

C3 

28.34 

34.19 

71.49 

C4 

50.39 

29.29 

-3.58 

c5 

-72.67* 

-86.77 

-35.32 

C6 

-29.22 

5.55 

-9.36 

C7 

21.93 

15.05 

9.18 

c8 

8.82 

0.78 

* 

4 -DAY  SAMPLE  BASE 


Band  1 

Band  2 

Band  3 

Sample  Size 

206 

227 

216 

Co 

9226.42 

8702.09 

7289.23 

cl 

15.01 

7.09 

* 

C2 

-36.75 

-20.98 

-16.37 

C3 

31.79 

49.50 

66.94 

C4 

52.29 

16.47 

* 

C5 

-77.84 

-83.39 

-33.79 

C6 

-27.24 

4.64 

-9.98 

C7 

26.96 

15.72 

5.07 

c8 

3.44 

0.41 

2.85 

* 

Indicates  an  insignificant  contribution  to  the  analysis  of 
variance  by  this  predictor. 
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The  FNWC  regression  coefficients  listed  in  Tables  5  and 
6  lend  support  to  the  observations  from  the  NESC  regression 
coefficients,  in  that  the  coefficients  from  the  FNWC 
dependent  test  vary  sufficiently  between  latitude  bands  to 
justify  using  separate  bands.   Furthermore,  the  differences 
between  the  three-day  samples  and  four -day  samples  are  again 
small;  and  the  coefficients  corresponding  to  Channels  3,  4, 
and  5  appear  to  be  weighted  more  heavily  than  the  others. 

Table  7  lists  the  values  for   r2   and  standard  errors 
for  each  of  the  FNWC  dependent  tests  listed  in  Tables  5  and 
6.   As  in  the  NESC  tests,  the  values  for   R   were  higher  in 
Band  3,  standard  errors  were  generally  larger  for  the  300  mb 
tests  than  for  the  500  mb  tests,  and  the  values  for   R^ 
were  higher  for  the  300  mb  tests  than  for  the  500  mb  tests. 
A  comparison  of  the  three-  and  four -day  sample  regressions 
for  both  the  300  mb  and  500  mb  cases  revealed  no  significant 
differences  in  corresponding  values  of  R^,  when  subjected  to 
to  a  t -test. 


34 


TABLE  7.   R2  and  Standard  Errors  (gpm)  for  FNWC 
Dependent  Regression  Equations  for 
500  mb  and  300  mb  Three-  and  Four -day 
Dependent  Samples. 


3 -Day  Sample 

R2 
Std.  error (m) 


500 

mb 

Rand  1 

Band  2 

Band  3 

.7910 

.7954 

.9014 

62.53 

84.55 

72.48 

4 -Day  Sample 


R' 


Std.  error (m) 


3 -Day  Sample 


Std.  error (m) 


.7787 

.7899 

63.79 

83.69 

300 

mb 

Band  1 

Band  2 

.8575 

.8361 

84.38 

105.99 

.8989 
73.21 


Band  3 

.9235 

84.32 


4 -Day  Sample 

R2 
Std.  error (m) 


.8459 
85.77 


.8328 
104.26 


.9222 
85.07 


35 


C.   COMPARISONS  BETWEEN  NESC  AND  FNWC  DEPENDENT  REGRESSION 
ANALYSES 

If  one  compares  the  results  from  Tables  2,  3,  5,  and  6, 
channels  3,  4,  and  5  appear  to  be  weighted  more  heavily  than 
the  other  channels  for  both  the  NESC  and  FNWC  dependent 
tests.   Table  8  combines  the  results  of  Tables  4  and  7  to 
facilitate  comparison  of  NESC  and  FNWC  dependent  regression 
test  results.   Here,  only  comparative  values  of  R2  by  the 
various  dependent  stratifications  are  considered.   If 
standard  errors  are  required,  reference  should  be  made  to 
Tables  4  and  7. 

Band  3  gave  the  highest  values  of  R  ,  averaging  8.7% 
higher  over  all  dependent  regression  systems  than  Bands  1 
and  2,  varying  from  5.6%  to  12.0%  higher  for  individual 
comparisons.   The  300  mb  values  of  R2   averaged  4.65% 
higher  overall  than  the  500  mb  values  of  R2.  with  the  NESC 
300  mb  values  of  R2   averaging  5.01%  higher  than  the  NESC 
500  mb  values,  and  the  FNWC  300  mb  values  of  R2   averaging 
4.38%  higher  than  the  FNWC  500  mb  values. 

Although  quite  close  in  most  instances,  in  all  cases  the 
NESC  values  of  R2   were  slightly  higher  than  the  correspond- 
ing values  of  R2   from  the  FNWC  tests,  averaging  2.01% 
higher  for  the  NESC  tests.   A  closer  look  reveals  an  average 
of  1.69%  smaller  values,  or  "shrinkage,"  for  the  FNWC  500  mb 
tests  as  compared  to  the  corresponding  NESC  tests.   An 
average  shrinkage  of  2.32%  occurred  for  the  300  mb  FNWC 
'tests  as  compared  to  the  NESC  tests.   The  individual 
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shrinkages  of  R2  for  each  band,  averaged  over  both  the  three- 
and  four-day  data-samples  is  shown  for  the  FNWC  tests  rela- 
tive to  the  NESC  corresponding  tests: 

Band  1     Band  2     Band  3 
500  mb       1.96%      2.61%      0.49% 

300  irib      2.60%      3.13%      1.21% 

2 

Bands  3  shows  the  smallest  shrinkage  of  R   values  for  FNWC  as 

compared  to  the  NESC  dependent  tests.   That  the  NESC  values 
of  R2  are  higher  than  the  FNWC  values  is  not  too  surprising 
since  the  NESC  data  was  originally  the  result  of  essentially 
an  eight-predictor  equation,  the  precise  form  of  which  was 
not  known  in  this  study. 

A  comparison  of  three-  and  four -day  overall  regression 
systems  showed  that  the  three-day  values  of  R  were  .46% 
higher  than  the  four-day  values.   Since  this  is  such  a  small 
difference,  it  is  likely  to  be  caused  only  by  the  additional 
data  samples  used  in  converting  the  particular  three-day 
sample  to  a  four -day  sample,  and  the  difference  was  not 
considered  significant. 

Thus  far,  it  is  seen  from  the  dependent-data  tests  on 
NESC  and  FNWC  samples  that  the  results  are  generally  good, 
considering  the  r2  values  averaged  .8579,  with  averages  of 
.8679  for  NESC  and  .847  8  for  FNWC.   From  the  paragraphs 
preceding  this  one,  it  appears  that  the  dependent  tests 
produced  better  results  in  Band  3  than  the  other  bands, 
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better  results  at  300  nib  than  at  500  nib  and  better  results 
for  NESC  than  for  FNWC  specification. 

It  should  be  pointed  out  that  the  NESC  and  FNWC  height 
data  were  highly  correlated.   The  500  mb  values  of  R^  ranged 
from  .9267  to  .9699  and  the  standard  errors  varied  from  35.26 
meters  to  43.15  meters.   The  values  of  R^  for  the  300  mb  case 
ranged  from  .9244  to  .9707  and  the  standard  errors  varied 
from  51.46  meters  to  59.64  meters. 
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IV.   INDEPENDENT  SAMPLE  TESTS 

A.   VERIFICATION  PROCEDURE  FOR  THE  INDEPENDENT  DATA  TESTS 

Each  of  the  dependent  regression  equations,  stratified  by 
NESC  and  FNWC,  by  latitude  bands,  and  for  the  different  sizes 
of  data-bases  in  Chapter  III,  was  tested  using  independent 
data  samples.   The  equations  developed  from  the  three-day 
samples  were  tested  with  independent  data  from  a  fourth  day, 
for  times  centered  at  12  hours  and  24  hours  later.   The 
regression  equations  developed  from  the  four -day  samples  were 
likewise  tested  for  times  centered  12  and  24  hours  later 
using  the  equations  developed  from  data  collected  over  24-27 
December,  1970.   The  corresponding  four -day  independent  tests 
utilized  data  from  the  28th  (28/00  GMT  and  28/12  GMT,  i.e., 
12  and  24  hours  later) .   In  addition,  the  12-  and  24-hour 
tests  provided  homogeneous  statistics  and  were  pooled  to 
provide  a  larger  number  of  independent  data  samples  for  each 
band  and  height. 

The  results  of  the  independent  tests  generate  equations 
of  the  general  form 

ZN(5,3)  =  DQ  +  0-^(5,3)  (11) 

and  ZF(5,3)  =  FQ  +  F-jZpfS^)  (12) 

where   ZN(5,3)   and   Zp(5,3)   represent  the  height-estimators 
from  the  regression  equation  (Equation  8) ,  DQ   and   FQ   are 

regression  constants,  and   D-j_   and   F-^   the  regression 
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coefficients.   ZF(5,3)   are  the  observed  independent-test 
heights  determined  by  interpolation  from  the  data  tapes  of 
FNWC  for  the  500  mb  and  300  mb  levels.   All  independent 
Zn(5,3)  were  the  result  of  application  of  the  NESC  update 
regression  to  the  independent  set  of  scan-spot  corrected 
radiances.   Regression  analyses  for  these  tests  were  similar 
to  the  regression  techniques  applied  in  Chapter  III,  in  that 
the  newly  computed  height  Z   is  now  being  compared  by 
regression  methods  to  the  known  height  at  each  scan-sounding 
location.   The  independent' test  results  are  presented  in  a 
slightly  different  manner  in  Tables  9,  10,  11,  and  12.   The 
mean  and  standard  deviation  for  the  NESC  and  FNWC  verifica- 
tions, and  the  values  for  R^  and  standard  error  resulting 
from   ZN(5) ,  ZN(3),  ZF(5),  and  ZF(3)   are  listed  in  Tables  9, 
10,  11,  and  12  for  the  three-  and  four -day  independent  tests. 


41 


TABLE  9.   Verification  Resulting  From  the  Use  of 

Equation  (8)  Upon  500  mb  NESC  Independent 
Data. 


3 -DAY  TEST 


Band  1 


Band  2 


Band  3 


Band  1 


Band  2 


Band  3 


Sample 
Size 

Mean 

Std. 
Dev. 

R2 

Std. 
Error 

12 

hr 

23 

5753.26 

99.04 

.7983 

45.52 

24 

hr 

35 

5738.60 

121.72 

.7  867 

57.07 

12+24 

hr 

58 

5744.41 

112.59 

.7852 

52.64 

12 

hr 

31 

5464.45 

100.7  6 

.5428 

69.30 

24 

hr 

37 

5477.38 

154.80 

.7652 

76.08 

12+24 

hr 

68 

5471.48 

132.15 

.6934 

73.62 

12 

hr 

32 

5283.97 

215.61 

.8781 

76.54 

24 

hr 

39 

5293.51 

227.12 

.9015 

72.24 

12+24 

hr 

71 

5289.21 

220.49 

.8859 

75.01 

4 -DAY  TEST 


hr 

Sample 
Size 

28 

Mean 
57  60.27 

Std. 
Dev. 

R2 

Std. 
Error 

12 

90.95 

.7537 

46.00 

24 

hr 

57 

5739.58 

121.90 

.8308 

50.59 

12+24 

hr 

85 

5746.36 

112.52 

.8151 

48.67 

12 

hr 

33 

5513.99 

140.14 

.7093 

76.77 

24 

hr 

50 

5432.87 

164.26 

.7499 

83.00 

12+24 

hr 

83 

5465.11 

159.32 

.7509 

80.00 

12 

hr 

26 

5290.45 

202.12 

.7685 

99.26 

24 

hr 

47 

5286.16 

172.66 

.8281 

72.38 

12+24 

hr 

73 

5287.69 

182.30 

.7959 

82.94 
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TABLE  10.   Verification  Resulting  from  the  Use  of 

Equation  (8)  Upon  300  mb  NESC  Independent 
Data. 


3 -DAY  TEST 


Sample 
Size 

Mean 

Std. 
Dev. 

R2 

Std. 
Error 

Band 

1 

12 

hr 

23 

9427.13 

166.40 

.8980 

54.39 

24 

hr 

35 

9398.46 

190.90 

.8526 

74.40 

12+24 

hr 

58 

9409.82 

180.62 

.8658 

66.76 

Band 

2 

12 

hr 

31  - 

8941.29 

157.58 

.7172 

85.23 

24 

hr 

37 

8981.95 

202.08 

.8027 

91.04 

12+24 

hr 

68 

8963.41 

182.96 

.7704 

88.34 

Band 

3 

12 

hr 

32 

8682.44 

293.24 

.9084 

90.23 

24 

hr 

39 

8695.46 

303.44 

.9334 

79.33 

12+24 

hr 

71 

8689.59 

296.83 

.9207 

84.19 

4 -DAY  TEST 

1 

12 

hr 

Sample 
Size 

28 

Mean 
9418.45 

Std. 
Dev. 

160.38 

R2 

Std. 

Error 

Band 

.8348 

66.42 

24 

hr 

57 

9409.62 

203.85 

.8718 

7  3.64 

12+24 

hr 

85 

9412.46 

189.70 

.8615 

71.03 

Band 

2 

12 

hr 

33 

9012.87 

206.88 

.7994 

94.15 

24 

hr 

50 

8923.04 

234.25 

.8137 

102.15 

12+24 

hr 

83 

8958.68 

226.81 

.8148 

98.22 

Band 

3 

12 

hr 

26 

8696.30 

266.23 

.8301 

111.99 

24 

hr 

47 

8696.93 

242.90 

.8827 

84.13 

12+24 

hr 

73 

8696.62 

249.60 

.8603 

93.94 
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TABLE  11.   Verification  Resulting  From  the  Use  of 

Equation  (8)  Upon  500  FNWC  Independent 
Data. 


3 -DAY  TEST 


Band  1 


Band  2 


Band  3 


Band  1 


Band  2 


Band  3 


hr 

Sample 
Size 

23 

Mean 

Std. 
Dev. 

91.62 

R2 

Std. 

Error 

12 

5736.63 

.7572 

46.21 

24 

hr 

35 

5724.33 

130.80 

.7103 

71.46 

12+24 

hr 

58 

5729.20 

116.11 

.7130 

62.75 

12 

hr 

31 

5468.17 

123.09 

.6163 

77.56 

24 

hr 

37 

5473.50 

159.20 

.7401 

82.31 

12+24 

hr 

68 

5471.07 

142.86 

.6950 

79.49 

12 

hr 

32 

5275.86 

208.00 

.8440 

83.52 

24 

hr 

39 

5291.72 

229.58 

.8863 

78.45 

12+24 

hr 

71 

5284.57 

218.72 

.8627 • 

81.62 

4 -DAY  TEST 


hr 

Sample 
Size 

28 

Mean 
57  50.87 

Std. 
Dev. 

R2 

Std. 

Error 

12 

80.03 

.8329 

33.34 

24 

hr 

57 

5728.56 

120.04 

.8054 

53.43 

12+24 

hr 

85 

5735.90 

108.52 

.8104 

47.53 

12 

hr 

33 

5522.73 

144.47 

.6730 

83.94 

24 

hr 

50 

5427.71 

171.16 

.7143 

92.43 

12+24 

hr 

83 

5465.49 

166.85 

.7170 

89.31 

12 

hr 

26 

5280.39 

198.52 

.6612 

117.93 

24 

hr 

47 

5274.83 

180.50 

.7782 

85.95 

12+24 

hr 

73 

5276.81 

185.76 

.7297 

97.26 
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TABLE  12.   Verification  Resulting  From  the  Use  of 

Equation  (8)  Upon  300  mb  FNWC  Independent 
Data. 


Band  1 


Band  2 


Band  3 


Band  1 


Band  2 


Band  3 


3- 

-DAY  TEST 

hr 

Sample 
Size 

23 

Mean 
9396.04 

Std. 
Dev. 

160.42 

R2 

Std. 
Error 

12 

.8465 

64.34 

24 

hr 

35 

9378.17 

195.55 

.7  643 

96.36 

12+24 

hr 

58 

9385.25 

181.16 

.7883 

84.10 

12 

hr 

31  . 

8934.96 

187.09 

.7281 

99.23 

24 

hr 

37 

8975.52 

217.55 

.7903 

101.03 

12+24 

hr 

68 

8957.03 

203.76 

.7673 

99.04 

12 

hr 

32 

8665.11 

282.86 

.8915 

94.70 

24 

hr 

39 

8694.04 

308.89 

.9157 

90.88 

12+24 

hr 

71 

4- 

8681.00 
-DAY  TEST 

295.70 

.9022 

93.13 

hr 

Sample 
Size 

28 

Mean 
9410.61 

Std. 
Dev. 

144.12 

R2 

Std. 
Error 

12 

.7820 

68.57 

24 

hr 

57 

9392.98 

200.90 

.8310 

83.34 

12+24 

hr 

85 

9398.79 

183.45 

.8198 

78.34 

12 

hr 

33 

9005.81 

217.27 

.7295 

114.81 

24 

hr 

50 

8913.67 

244.64 

.7718 

118.07 

12+24 

hr 

83 

8950.30 

237.16 

.7  641 

115.91 

12 

hr 

26 

8678.84 

267.39 

.7858 

126.31 

24 

hr 

47 

8679.63 

257.69 

.8367 

104.31 

12+24 

hr 

73 

8679.34 

259.33 

.8188 

111.18 
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B.   SUMMARY  OF  INDEPENDENT  TEST  RESULTS 

1.   Results  Internal  to  the  Independent  Test 
Table  13  lists  the  combined  12  and  24  hour 
independent-test  values  of  R^   abstracted  from  Tables  9, 
10,  11,  and  12  to  facilitate  comparisons  of  test-methods 
contained  in  those  tables.   Again  our  emphasis  is  on  R2   as 
a  measure  of  the  efficiency  of  the  verification  since  this 
statistic  has  the  effect  of  normalizing  the  variable  standard 
deviations  in  Tables  9,  10,  11,  and  12. 

Band  3,  as  in  the  dependent  tests,  gave  the  highest 
values  of  R  ,  averaging  6.99%  higher  than  the  other  bands. 
The  300  mb  values  of  R  ,  again,  as  in  the  dependent  cases, 
were  higher  than  the  500  mb  values,  averaging  5.83%  higher. 
The  NESC  values  of  R2 ,  as  compared  to  the  FNWC  values  of  R2, 
were  higher,  as  before,  by  an  average  of  3.60%  over  all  of 
the  independent  tests.   However,  the  amount  of  shrinkage 
between  the  NESC  and  FNWC  within  the  independent  tests  was  a 
minimum  in  Band  2,  as  compared  to  the  minimum  values  of 
shrinkage  in  Band  3  for  dependent  tests.   This  phenomenom  is 
due  to  the  fact  that  the  independent  test  of   RN    fell  off 
considerably  in  Band  2,  whereas  the  corresponding  shrinkage 

was  somewhat  less  with  FNWC  independent  tests  in  Band  2.   The 

2 

three-day  estimator  tests  produced  average  values  of  R 

that  were  .7  6%  higher  than  the  corresponding  four -day  values 
for   R  .   This,  again,  as  in  the  dependent  tests,  is  a  small 
difference  and  is  considered  to  be  due  to  the  small  sample- 
size,  or  to  the  particular  data-sample  added  to  obtain  the 
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four -day  sample,  and  therefore  of  little  statistical 
significance. 

2.   Comparison  of  Independent-  With  Dependent-Data 
Specification. 

From  the  independent  tests,  the  same  basic 
conclusions  as  for  the  dependent  tests  are  to  be  drawn  since 
the  results  are  quite  similar.   The  independent-test  values 
of  R2   were  generally  good  with  an  average  of  .8004,  an 
average  of  .8184  for   R2   from  the  independent  NESC  tests 
and  an  average  of  .7824  for  the  FNWC  tests.   The  highest 
values  of  R2  were  in  Band  3,  highest  at  300  mb,  and 
highest  for  the  NESC  tests. 

In  proceeding  from  the  dependent  tests  to  the 
independent  tests,  some  loss  of  accuracy  is  to  be  expected. 
Values  of  R2   again  will  be  used  for  comparisons.   Table 
14  is  derived  from  Tables  8  and  13  to  present  an  easier 
comparison  of  the  corresponding  values  of  R2   from  the 
dependent  tests  and  the  independent  tests,  with  the  value 
of  shrinkage  indicated  also.   In  most  cases  the  shrinkage 
of  the   R2   values  from  the  FNWC  dependent  tests  to  the 
independent  tests  was  slightly  larger  than  the  shrinkage 
for  the  similar  NESC  tests.   The  FNWC  shrinkage  averaged 
1.59%  more  than  the  NESC  shrinkage. 

The  first  group  of  numbers  in  Table  14  are  averaged 
over  all  bands  and  over  both  three-  and  four-day  coeffi- 
cient estimators.   In. this  group,  the  values  of   R2   for 
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the  NESC  300  mb  tests  show  the  lowest  relative  shrinkage  of 
the  four  shown  in  the  group.   This  might  have  been  expected 
since  the  300  mb  NESC  test  results  from  both  the  dependent 
and  independent  tests  showed  better  results  than  the  others. 

The  second  group  of  numbers  in  Table  14  is  averaged 
over  the  500  mb  and  300  mb  level  independent  tests,  and 
over  three-  and  four-day  coefficient  estimators   Z.   This 
second  group  shows  relatively  small  shrinkages  in  Band  1. 
This  could  possibly  be  due  to  the  fact  that  in  low  latitudes, 
the  radiance  correction  p'rocedure  is  somewhat  more  consis- 
tent than  in  mid-latitudes,  and  may  be  affecting  the  amount 
of  shrinkage  in  a  favorable  manner. 

The  third  group  of  numbers  in  Table  14  is  averaged 
over  all  three  latitude  bands,  and  shows  the  relative 
shrinkages  by  use  of  the  Equation  (8)  based  upon  the  three- 
and  four -day  dependent  coefficient  matrices.   The  amount  of 
shrinkage  of  the  three-day  independent  tests  over  the  four- 
day  tests,  or  vice-versa,  is  not  very  obvious,  even  though 
one  might  expect  that  the  four -day  tests  should  be  more 
representative  due  to  a  larger  input  of  data.   From  the 
findings  of  this  study,  it  has  consistently  been  determined 
that  tests  on  independent  data  using  the  four-day  data-base 
estimators  are  on  the  average  slightly  less  effective  in 
R^   than  when  the  corresponding  three-day  data  base  was 
used. 
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TABLE  14.   Summary  of  Dependent  to  Independent 
Test  Shrinkage. 


Test  Being 
Compared 

Dependent 
Value  of  R2 

Independent 
Value  of  R2 

Average 
Shrinkage 

ZN(5) 

.8428 

.7878 

5.50% 

ZF(5) 

.8259 

.7547 

7.12% 

ZnO) 

.8929 

.8489 

4.40% 

ZF(3) 

.8697 

.8101 

5.96% 
avg.  5.7  5% 

ZN  Band 

1 

.8411   * 

.8320 

.91% 

ZF  Band 

1 

.8183 

.7829 

3 .  54% 

ZN  Band 

2 

.8423 

.7572 

8.51% 

ZF  Band 

2 

.8136 

.7359 

7.77% 

ZN  Band 

3 

.9201 

.8657 

5 .  44% 

Zp  Band 

3 

..9116 

.8283 

8.33% 

avg .    5.75% 


Zn(5)  3-day 

ZN(5)  4-day 

ZN(3)  3-day 

ZN(3)  4-day 

ZF(5)  3-day 

ZF(5)  4-day 

ZF(3)  3-day 

ZF(3)  4-day 


8441 
8415 
8945 
8912 
8293 
8225 
87  24 
8670 


.7882 

5.59% 

.7873 

5.42% 

.8523 

4.22% 

.8455 

4.57% 

.7569 

7.24% 

.7524 

7.01% 

.8193 

5.31% 

.8009 

6.61% 

* 

avg.  5.75% 
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V.   CONCLUSIONS 

It  has  been  shown  in  this  study  that  regression  equations 
developed  from  relatively  small  data  bases  of  heights  and 
corrected  SIRS-B  radiances  can  be  used  to  compute  heights  on 
independent  samples  with  relatively  good  results.   It  is  felt 
that  further  testing  should  be  conducted  at  other  levels 
using  comparisons  of  other  samples  to  test  the  results 
obtained  here.   If  the  resources  are  available,  larger  data 
bases  should  be  tested  in  a  similar  manner.   Direct  applica- 
tion of  similar  procedures  could  be  of  possible  use  in  the 
update  of  FNWC  analyses  for  standard  levels  in  areas  of 
sparse  data,  or  even  aboard  ships  with  weather  offices  to 
provide  real-time  height  calculations  over  an  operational 
area.   In  the  latter  case  the  advantage  of  a  small  data  base 
lies  in  the  relatively  small  amount  of  core  storage  necessary 
for  the  type  of  computer  calculations  possible  in  shipboard 
operations. 

Although  the  NESC  height  data  provided  slightly  better 
regression  results  than  the  FNWC  height  data,  this  was  to  be 
expected  since  the  NESC  heights  were  computed  using  an  eight- 
predictor  regression  equation,  and  therefore  were  internally 
more  consistent.   The  results  obtained  by  using  FNWC  data 
for  dependent  and  independent  regression  analyses,  parallel- 
ing the  procedures  applied  to  the  NESC  data,  lends  support 
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to  the  overall  concept  of  developing  regression  equations  to 
compute  geopotential  heights  using  SIRS-B  corrected  radiances 

Improvements  of  the  results  should  be  expected  if 
longer-term  data  samples  are  used.   Improvements  should  also 
be  expected  if  the  width  of  the  latitude  bands  were  smaller, 
and  if  the  radiance  data  were  subjected  to  gross  error  checks 
before  entering  into  regression  analyses.   Also,  a  fact  to 
be  considered  is  the  tendency  of  the  SIRS-B  channel  response 
to  degrade  and  require  updated  calibration  after  being  in 
orbit  for  a  period  of  time".   If  degradation  errors  were 
linear  in  time,  the  regression-climatology  should  provide 
update  corrections,  but  the  normal  case  is  that  degradation 
occurs  non-uniformly  in  the  different  channels.   The  best 
scan-soundings  are  normally  obtained  soon  after  the  satellite 
is  launched  into  its  orbit  about  the  earth,  before  appreci- 
able system-degradation  can  occur. 
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APPENDIX 
GRIDDING  AND  INTERPOLATION  SCHEME 
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